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ABSTRACT 

In a power spectrum of the Sun's acoustic waves, the p modes have distinctly 
skewed frequency profiles. Furthermore, the asymmetry is observed to have the 
opposite sign in power spectra made from line-of-sight velocity and continuum 
intensity. The asymmetry and its reversal in sign has previously been explained 
using a combination of mechanisms that involve the acoustic source. A localized 
acoustic source within an acoustic cavity naturally generates asymmetric profiles 
through wave interference; however, the sign of the asymmetry due to this mech- 
anism is identical for all observables. The reversal of the asymmetry between 
velocity and intensity observations has been attributed to the visibility of the 
source itself (i.e., "correlated noise"). In this paper, I will show that asymmetry 
generated by a localized source can be interpreted as either a wave interference 
effect in physical space, or a mode interference effect in spectral space. I advo- 
cate a new mode-fitting procedure based on this new interpretation, whereby the 
complex phases of all the modes determine the mode asymmetries. Further, I 
suggest that information about the acoustic source function is encapsulated in 
the amplitude of each mode, and present a scheme by which the source func- 
tion can be obtained from measured mode amplitudes by standard helioseismic 
inversion techniques. 

Subject headings: Methods: data analysis — Sun: helioseismology — Sun: oscillations — 
Waves 



1. Introduction 



In power spectra of solar acoustic oscillations the frequency profiles of the p modes 
are skewed and markedly non-Lorentzian. This proper ty was first discove red using Ca II 
K observations made from the geographic South Pole (jPuvall et al.l Il993l ) and were later 
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confirmed (iNigam k Kosovichevlll998h with modem spacecraft-based telescopes using data 



from the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory 
(SOHO). In both sets of observations, a puzzling fact was revealed: for modes of the same 
frequency and harmonic degree, the asymmetry or skewness of the line profile has opposite 
signs for observations made in the continuum intensity and line-of-sight Doppler velocity. 
For both types of observations, the lowest-frequency modes are the most severely skewed 
and the asymmetry appears to be only a weak function of harmonic degree; however, for 
modes of high harmonic degree (/ > 150) this statement is difficult to verify because the 
mode peaks are highly blended with the spatial leaks from modes of nearby degree, forming 
a broad ridge of power. Figures [1] and [2] illustrate the asymmetry for both high and low 
degree. From the figures it is evident that the frequency profile of each individual mode is 
asymmetric and, at high degree, the blended ridge is also asymmetric. 

The physical mechanism for this skewness depends somewhat on whether we are con- 
sidering the asymmetry of each mode profile or the asymmetry of a blended ridge. In the 
former case, the asymmetry of each mode is now thought to arise from a combination of 
effects due to the nature and location of the acoustic wave source. At high harmonic degree, 
the asymmetry of the p mode ridge is the result of the skewness of each component mode 
profile as well as an uneven distribution of power amongst the spatial leaks. The wave source 
can introduce skewness in a variety of ways. The first mechanism that I discuss here relies 
in detail on the location of the acoustic source. The primary source of acoustic waves in the 
Sun is granular convection. Short-lived convective events of relatively high Mach number 
prolifically emit sound waves. Since granulation is confined to a narrow band just below the 
photosphere, the acoustic source is highly localized in depth. Waves generated in this layer 
propagate both upwards and downwards. Those that propagate downwards are eventually 
refracted back upwards at the bottom of the acoustic cavity by the sound speed gradient. 
When this initially downward propagating wave returns to the excitation layer, it can de- 
structively interfere with the wave that was initially upward propagating. At particular 
frequencies, complete destructive interference occurs and zero wave amplitude is observed 
above the excitation layer (which is where the spectral lines used in actual observations are 
formed). The frequencies at which the interference is complete are as ymmetrically placed 



between the mode ridges, thereby leading to asymmetric mode profiles (lGabriellll992l . Il995 



Roxburgh fc Vorontsov 1995 : Abrams fc Kumar 1996 : Rast fc Bogdan 19981 ). 



The second mechanism due to the source arises because the convective events that 
generate sound waves are themselves regions of intense temperature fiuctuation and high 
velocity. The sources therefore provide a contribution to the power in either intensity or 
velocity observations. Since the acoustic waves are correlated with the source that generates 
them, the visibility of the source can lead to power that varies significantly over the line 
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profile and cause skewness. This mechanism has been dubbed "correlated noise" and has 
been suggested as the reason why the asymmetry is opposite for intensity and v elocity spectra 
( iRoxburgh &: Vorontsovl 119971 : iNigam fc KosovichevI Il998l : iNigam et al.l Il998[ ) . The current 
understanding is that in order to reproduce the asymmetry of both intensity and velocity 
power spectra and to explain the rapid variation with frequ ency of the phase of the cross- 
spectra between intensity and velocity (lOliviero et al.l Il999l). correlated noise must provide 



significant power to both types of spectra (e.g., iNigam et al.lll999l : ISkartlien fc RastI 12000 
Jefferies et aLll2003h . 



Numerical simulations of granulation have provided a counter explanation (iGeorgobiani et al. 



20031 ) . In these simulations, the velocity and intensity fluctuations have power spectra with 
asymmetric mode profiles. However, if the two quantities are sampled at the same geometric 
height the sign of the asjTiimetries are identical. Only after artificial Dopplergrams and 
intensity images have been constructed using radiative transfer do the two types of power 
spectra exhibit opposite asymmetries. Therefore, these simulations suggest that details of the 
radiative transfer, such as the different heights at which the two observables are formed and 
the nonlinear interaction between wave fluctuations and the radiation field, are responsible 
for the opposite skewnesses. 

Unfortunately, many mode-fitting schemes used to r neasure p mode freq uencies assume 
that the line profile is symmetric, a Lorentzian in fact. iDuvall et al.l (119931 ) estimate that 
failing to account for the asymmetry in the mode's frequency profile can lead to mismeasure- 
ment of mode frequencies by as much as one part in 10^, a significant error for the purposes of 
helioseismic inversions. This is probably the source of the systematic offset (~ 0.1 /xHz) that 



was o bserved between frequencies obtained from intensity and velocity spectra (IToutain et al. 



19971 ). To account for such systematic shifts, researchers ha ve suggested a variety of asym 



metr i c model profiles with which to fit the observations fe.g.. lDuyall et al 



1993 



Nigam et al. 



19981 : iRosenthallEggsl ) . Due to its practicality, the formulation of lNigam et al.l (119981 ) is prob- 
ably the mo st widely used, having been adapted for both the meas urement of global mode 



frequencies (iToutain et al.l 



analysis spect ra (iBasu &: Antia 



1998 



Reiter et al 



1999 



2002 



Basu et al 



2001 



KorzennikI 120051) and the f itting of ring- 



Tripathy et al.l |2009| ) . In essence 



the profile of iNigam et al.l (119981 ) is a Lorentzian times a polynomial in frequency, where 
the polynomial has been calculated by performing a low order expansion about the mode 
frequency ojq, 
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In this formula, 7 is the hnewidth, A is the power amphtude, and s is an asymmetry param- 
eter. 

This expression is formally valid under fairly limited conditions. For example, it has 
been derived as a small argument expansion in terms of sx, which requires that asymmetry 
is weak and more restrictively that we confine our attention to the core of the line and avoid 
the wings. One clear sign of this restriction is that the profile isn't integrable, and does not 
contain a finite amount of energy. Second, it ignores the presence of nearby modes which 
may contribute to the power. This is a particularly onerous assumption for high degree 
modes for which the line widths are sufficiently large that the mode and the neighboring 
spatial leaks blend into a ridge. Figure [1] illustrates this later point quite well. Each ridge is 
a blending of mode peaks and the ridges are wide enough that they have begun to overlap. 

Our goal here is to generate a fitting function that takes the skewness of the p-mode 
frequency profiles into account while maintaining its validity over a wide range of conditions. 
I do so by using a modal expansion of the wavefield to represent the power spectra using 
a small number of free parameters. This expansion is accomplished in a general manner 
that doesn't depend on details of the atmosphere or the mechanisms for wave damping and 
excitation. In the next section. Section |2l I discuss the driving and damping of acoustic p 
mode oscillations in broad terms. In Section [3l I derive the modal expansion of the wave 
field and compute the resulting power spectrum. In Section |U I promulgate a fitting profile 
of the spectrum that treats the asymmetry as a mode interference phenomenon. Finally, in 
Section [5l I discuss my findings and suggests ways in which properties of the acoustic source 
could be deduced through inversion of the fitted complex mode amplitudes. 



To avoid unnecessary complexity, consider a neutrally-stable, plane-parallel atmosphere — 
a reasonable approximation for modes of intermediate degree and larger (/ > 20) which are 
trapped within the solar convection zone. Acoustic waves in such an atmosphere obey the 
following wave equation: 



where the sound speed c and acoustic cut-off frequency Uc are functions of height z. Tech- 
nically, the damping rate 7 may also vary with height, but for simplicity I will restrict my 
attention to a constant damping rate. The quantity ip can represent any wave variable (with 



2. 



Driven Acoustic Oscillations 
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a suitable definition for tlie acoustic cut-off frequency); liowever, in order to relate our find- 
ings to observed power spectra made from Dopplergrams, I will treat ip as proportional to the 
vertical component of the velocity. The function S{x,t) represents the density of acoustic 
sources. 



Since the atmosphere is invariant in the horizontal direction and in time, Equation (12.11) 
can be Fourier transformed trivially to obtain an ODE in only the height coordinate z, 



+ ^^-V{z; k) 



dz'^ ' c^{z) 
where I adopt several definitions, 



■ip{z; k, u) = S{z; k, u) 



(2.2) 



Viz; k) 

and use the following Fourier convention: 
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(2.5) 
(2.6) 



The frequency fl is inherently complex and the acoustic potential V is purely real. Note, 
that if the complex wave frequency u = Ur + ioJi has an imaginary part Ui = —7, then the 
frequency is purely real with = ^4;^"+^. This allows the solution to the damped 
problem to be generated by solving an equivalent undamped problem and then shifting the 
resulting frequency, in both its real and imaginary components, by an amount dependent on 
the damping rate, u = 1/^^—7^ — i'~f. 



2.1. Green's Function 



I commence by solving for the spectral Green's function G. For a point source located 
at the height z', the Green's function satisfies 
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G{z,z';k,u) = 5{z-z') . (2.7) 

The Green's function can be found in the standard way. If I denote ip- as the homogeneous 
solution that obeys the boundary conditions deep in the atmosphere (as z — )■ — oo) and 
similarly write ijj^ as the homogeneous solution that obeys the boundary condition high in 
the atmosphere (as z — )■ +cxd), the Green's function may be written in the following form: 

^/ / , N "^(z, z': k,u) , , 

where W is the Wronskian of the two homogeneous solutions and is their product, 



— + — -V 



d^p+ dip_ 
dz dz 



+ 5 



ip^^z; k,ijj) ip^{z'; k,oj) if z < z' 
ip-{z'; k,ijj) iIj+{z; k,ijj) if z > z' 



(2.9) 
(2.10) 



Equation (12. 8 p has been used extensively in the p ast to explain the asymmetry of the p- 
mode frequency profiles in the Sun's acoustic spectrum (lGabriellll992l : [Roxburgh fc Vorontsov 
19951 : lAbrams fc Kumar! Il996l : iRast fc BogdanI Il998h . Since the driving of acoustic waves is 
thought to arise from granulation, which is localized in a thin layer just below the photo- 
sphere, I make the assumption that the height dependence of the source function is in fact 
just a delta function. 



S{z;k,uj)=6{z-z')Soik,u) . (2.11) 

This means that the Green's function provides the correct vertical behavior of the full wave- 
field. 

If one further assumes that the spectral line used to observe the wavefield forms at a 
height that overlies the driving layer {z > z'), the solution becomes 

^{z; k, u) = So{k, u) — ^ . 2.12 

The power spectrum of such a wavefield, observed at a height z in the atmosphere, is of 
course 
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i2 ir./,_ , .^,2\'^-{z';k,u)f \ilj+{z;k,u)f 



which may be thought of as having four separate contributions: (1) the strength of the source 
|S'o|^, (2) a contribution from the Wronskian |W|~^ that represents the modal structure of 
the acoustic cavity, (3) a contribution from the driving efficiency \iIj-{z'; k,u)f , and (4) a 
contribution from the height of observation {ip+iz] k, uj)f. The Wronskian possesses a zero at 
each mode frequency (which in the present formulation will be a complex frequency because 
of damping). These zeros generate Lorentzian peaks in the power spectrum at each mode 
frequency and by themselves these peaks are symmetric functions of frequency. The observed 
asymmetry arises from the driving efficiency which expresses how well the driver couples to 
any given mode. The driving efficiency possesses zeros wherever the homogeneous solution 
ip- has a node at the driving height. This occurs whenever the downward propagating wave 
generated by the source rebounds off the lower turning point of the acoustic cavity and 
returns to the driving layer with the exact phase needed to destructively interfere with the 
upward propagating wave also excited by the source. When this occurs, the solution vanishes 
for all heights above the driving layer. The zeroes in the driving efficiency fall in between 
the mode frequencies and are not distributed symmetrically about the mo de peaks; hence 



they generate asymmetry in the total power profile (for complete details see lRast fc Bogdan 



1998|). 



2.2. A Simple Example 

To illustrate the salient features of this behavior without undue complications, consider 
a simple "guitar-string" model. Waves are trapped between two heights in the atmosphere, 
z = and z = L, ai which the solution ip must vanish. Further, within this cavity let the 
sound speed be constant and let the acoustic potential vanish. The solutions are therefore 
purely sinusoidal with wavenumber K = Q/c, 



ij-{z;u) = sm{Kz) , (2.14) 
Mz;uj) = sm[K{L-z)] , (2.15) 
W(a;) = Ksm{KL). (2.16) 

If 1 assume that the driving height lies below the observation height, and further assume that 
the spectral dependence of the driver is white (i.e., 5*0 is a constant function of frequency), 
then the power spectrum of the solution is as follows: 
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P{co) = \S^ 







^sin^ {Kz') sin^ [K{L - z)] 
K'^ sin^ (KL) 



(2.17) 



The frequency dependence of this expression is illustrated in Figure [3] with the thick 
black curve. The individual contributions from the Wronskian, |i^'sin(i^'L)|~^ (red curve), 
the driving efficiency, sin^^Kz') (blue curve) and the observation height, sin^ [K{L — z)] 
(green curve), are also indicated. The line asymmetry arises because the spacing between 
the zeros in the driving efficiency is smaller than the spacing between modes (i.e., the spacing 
between the zeros of the Wronskian). The most asymmetric modes (those at low frequency) 
are those for which a zero in the driving efficiency almost coincides with a mode frequency. 
For the model presented in the figure, the height of observation is very close to the upper 
boundary of the acoustic cavity (as is true for helioseismic observations of solar p modes). 
This leads to a slow frequency variation in the power due to the height of formation; therefore, 
this term only provides a slow variation in the overall amplitude instead of skewness. 



While Equation (12.131) is useful for understanding why the Sun's p-mode profiles are 
observed to be asymmetric functions of frequency, this formulation for the solution ip and 
the resulting power spectrum P does not lend itself readily to the problem of mode fitting, 
whereby the mode parameters (i.e., frequency, lifetime and amplitude) are measured from an 
observed spectrum. Explicit use of Equation (12.131) requires prior knowledge of the homoge- 
neous solutions ip- and which pre-supposes a detailed understanding of the structure of 
the solar atmosphere (obtaining such an understanding may be the reason that one is mea- 
suring the mode frequencies in the first place). A more natural representation of the solution 
would be a modal decomposition, where the mode parameters are represented explicitly. 

If one ignores the existence of a global cut-off frequency and therefore neglects the 
possibility of leaky modes and continuous spectra, one may use the orthogonality properties 
of the eigenfunctions of the homogeneous problem. 



to perform a standard eigenfunction expansion of the solution, generating the following 
expression for the wavefield: 



3. Modal Decomposition 



oo 





— oo 
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n=l 



ipniz'; k) tpnjz; k) 

- niik) 



(3.2) 



When the atmosphere possesses a global cut-off frequency (as is true for the Sun's p modes), 
an eigenfunction expansion is still possible; however, the expressions are more complicated 
due to the existence of continuous spectra. To avoid obscufating the important issues with 
unneeded mathematics, I will consider such possibilities separately in Appendix |Al In Equa- 
tion (13. 2p the eigenvalues fl"^ and eigenfunctions ipn are purely real, whereas fi^ is complex. 
This formulation is completely equivalent to that shown previously in Equation f l2.12p : how- 
ever, in this treatment the contribution from each mode to the power is explicitly represented. 



P(k,u) 



\So{k,u) 



E 

n=l 



ipn{z'; k)iJn{z; k) 



(3.3) 



Note that each eigenfunction tpn is associated with two separate complex eigenfrequencies 
u±n{k) = ±Unik) — ^7, one with a positive real part (which I will denote with the positive 
index n) and one with a negative real part (denoted with the negative index —n). Therefore, 
the contribution from each eigenmode can be decomposed into its two frequency components. 



'ilj{z;k,uj) = So{k,uj) 



P(k,uj) 
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oo 
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\Sn(k,uj) 
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20Jn{k) UJ — (jJnik) ' 

i)n{z'\k)i)n{z;k) l-5no 



2uJn{k) 00 — OJn{k) 



(3.4) 
(3.5) 

(3.6) 



The summation over n now includes the negative indices, the variable a;„ (= — W-n) is the 
real part of the complex mode frequency, and the factor with the Kronecker delta 5^0 has 
been included to remove the n = term from the summation. 

The familiar representation of the power spectrum as the sum of modes, each with a 
Lorentzian profile, can be recovered by assuming that all of the modes are well-isolated 
from each other spectrally (the damping rate is much less than the spacing between modes). 
Therefore, mode interference cross-terms can be ignored and the power is just the incoherent 
sum of the power within each mode. 



- 10 - 



Pinc{k,u) = \So{k,u;)\ -— rry— ^ • (3-7) 

Of course, each Lorentzian profile in Equation fl3.7p is symmetric about the central mode 
frequency ujn{k)- The asymmetry that is observed in p-mode hne profiles comes from the 
mode interference terms that were neglected in Equation (13. 71) . 



3.1. The Simple Example Revisited 

The fact that interference between modes can be construed as the source of asymmetry 
in p-mode power profiles is aptly illustrated by the simple guitar-string model previously 
explored in Section 12.21 For this simple problem, the eigenvalues and eigenfunctions are 
defined by 





= KnC = 


nnc 
L ' 




= ±OJn - 


i'-f = z 
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(3.8) 
(3.9) 

(3.10) 



and the modal expansion is just a Fourier sine series. 



ij{z;uj) 



Pioo) 



L 



Esinimrz' /L) sm(nTcz/L) - 



n 



L2 



N sm(nn z' / L) sm(n7i z / L) - 



(3.11) 
(3.12) 



The black curve in Figure [3] represents the power from either representation, equations fl2.17p 
or f l3.12p . However, from the perspective of an eigenfunction expansion, the asymmetry in the 
frequency profile of each mode arises because of mode interference. Zeroes in the power are 
located at those frequencies for which complete destructive interference has occurred, where 
the interference in this interpretation is in spectral space between modes with overlapping 
power not in physical space between different wave components generated by the source. 
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Direct confirmation of tliis comes from plotting tfie incolierent sum of the power in all the 
modes, 

Pincioo) = T sm\nnz'/L) sin^nnz/L) ' , , . (3.13) 

In FigureOthis sum is shown as the orange curve and it can clearly be seen that the frequency 
profiles are symmetric, albeit modified by a slowly varying background (which is formally 
composed of the wings of many modes). 



4. Mode Fitting 

In a traditional procedure, to measure properties of a p mode such as its frequency, 
amplitude and lifetime, the power within a narrow band of frequencies centered on the 
target mode is fit with a model function. The range of frequency is selected with the goal 
of mitigating the contamination from other modes. However, since spatial leakage produces 
spectral peaks that cannot be fully isolated from the target mode (see Figured]), the fitting 
function is usually comprised of an incoherent sum of power peaks based on the properties of 
Equation f l3.7p . For the jth mode (where j indicates a particular combination of harmonic 
degree, azimuthal order and mode order) the fitting function might have the form, 

Pmodel(w) = B{UJ) + Vj{uj) + ^ Vj>{uj) (4.1) 

where each profile 'Pj{u) is a Lorentzian function with three parameters. 



In these expressions, i? is a background function that represents incoherent solar noise, 
and the three parameters that describe the Lorentzian are its mode frequency Uj, power 
amplitude Aj and line width Tj. The leakage is described by the leakage matrix Cjj' and 
the summation is performed over all modes that have a frequency within the fitting window 
and that have significant leakage. All of the free parameters in this formulation are treated 
as real quantities. 



Modern fitting procedures (e.g.. iKorzennikI |2005| ) often attempt t o account for t he ob - 
served skewness in the mode profiles by using the asymmetric profile of iNigam et al.l (119981 ) . 
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Equation (11. ip . Using the same notation presented here, Equation (ll.ip becomes 





( 




(4.3) 



which adds one additional parameter, the asymmetry parameter Sj, to each fitted peak. 

As mentioned previously, this commonly used skewed Lorentzian function is only a 
valid represention within the core of the line and it does not include the effects of mode 
interference (in spectral space). Further, the energy contained by any given mode (and 
hence the amplitude parameter Aj) is poorly defined since each mode's power profile lacks 
a finite integral. A globally valid expression that directly accounts for interference effects is 
suggested by the functional form of Equation f l3.3p . 



where the parameters all have the same meaning except that the mode amplitude aj and the 
leakage matrix Cjji are now complex quantities. Unlike classical fitting functions, this model 
function properly captures asymmetric mode profiles in both the wing and core because the 
effects of mode interference are present through the retention of the cross-terms between 
modes. The cost of accounting for the asymmetry is the need to fit one additional parameter 
per mode, the complex phase of the amplitude = \aj\e'^'^K 



I have suggested an alternate formula — Equation (14. 4p — for fitting p-mode frequency 
profiles that is based on an eigenfunction expansion of the Green's function. This formulation 
takes into account the skewness of the line profiles that is inherently produced by a localized 
acoustic source by including the effects of mode interference in spectral space. In addition, 
the fitting formula has the salutory features that it is valid for all frequencies across the mode 
profile (not just in the core of a line) and each mode contains a finite amount of energy. 

The suggested scheme can be applied to the fitting of global modes of relatively low 
harmonic degree (/ < 200) directly as envisioned. The asymmetry in a mode's profile is 
primarily the result of interference with nearby spatial leaks, which are most likely from the 
same mode order or ridge (except at very low harmonic degree / < 10). On the other hand. 




2 



(4.4) 



5. 



Discussion 
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the proposed fitting model is not immediately relevant to the fitting of ring-analysis spectra. 
The asymmetry observed in the p-mode ridges at large harmonic degree (Z > 200) are not 
produced by the same mechanisms that generate the asymmetry for global modes. Instead, 
the asymmetry is largely a result of the spatial window function (the ring-analysis version 
of the blending of a mode and the nearby leaks). 

Ring-analysis spectra are manufactured by Fourier transforming the wavefield observed 
on a small patch on the surface of the Sun. The small spatial domain translates into a 
broad spread of power in horizontal wavenumber. This widens the ridges and produces 
an effective line width that is much larger than one would expect based just on the mode 
lifetime. Furthermore, the resulting line profile depends on the window function applied 
during the Fourier analysis, the slope of the p-mode ridge in the dispersion diagram (dun/dk), 
and the variation in power along the ridge. Unlike in the global mode analyses where the 
effects of leakage have been considered in detail for over a decade, leakage has not been 
treated self-consistently in any ring-analysis studies. While the fitting formula presented here 
could be used to model the resulting ring-analysis spectra quite well with the asymmetry 
being produced by the interference between the ridges themselves, the underlying physical 
mechanism would be wrong. I leave it to a subsequent paper to include the effects of the 
leakage function on mode profiles. 



One novel property of the fitting function given by Equation (14. 4p is that the effect of the 
source is directly represented by the complex mode amplitudes. Therefore, the additional 
parameters have a physical significance that can be used to ferret out information about 
acoustic sources and the convection that engenders them. This can be seen most easily if I 
relax the assumption that the source is proportional to a delta function in height. If instead 
the source is distributed, as usual the wavefield is obtained through a convolution of the 
Green's function with the source. 



5.1. Measuring the Acoustic Source Function 





oo 




(5.2) 



Assuming that the source function is isotropic with a frequency variation that is slow com- 
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pared to the damping rate, the complex mode amphtude is therefore given by an integral 
over the source, 



anik) = / dz' ICn{z';k) S{z';k,Un) , (5.3) 



2uJn[k) 

where the set of /C„ are sensitivity kernels. Given a set of measured amplitudes, the coupled 
integral equations (15. 3 p could be inverted using standard helioseismic inversion procedures 
to obtain the source as a function of height in the atmosphere. 



5.2. Effects of Correlated Noise 



The theory and procedures that have been developed in this paper have so far ignored 
the existence of correlated noise. Previous studies have shown that correlated noise is a 
necessary ingredient to explain th e reversal of asymmetry betw' e en Doppler-velocity spectra 



and continuum-int ensity spectra (iRoxburgh fc Vorontsovl 119971 : iNigam fc KosovichevI Il998 



Nigam et al. 1998[). as well as th e rapid phase variation in velocity-intensity phase difference 



spectra ( 



Skartlien fc Rastll200Cll ). The modal ex pansion presented here c an easily be adapted 



Skartlien fc RastI 120001 ): one simply adds a 



to include the effects of correlated noise (e.g., 
complex correlated noise component, V'corr, to the wavefield before computing the power. 



-Pmodcl(c^) = Biuj) 



U — (jjj + 



(5.5) 



Assuming that the correlated noise is due to the visibility of the source itself, one would 
expect that the phase and amplitude of the correlated noise T/'corr would be a slowly varying 
function of frequency. Therefore, if one is fitting the power in a narrow frequency band, it 
can safely be assumed that the noise is a complex constant. Thus, in any fit to data the 
correlated noise adds only two additional parameters that must be fit. If one were fitting a 
broad range of frequencies, it might be necessary to parameterize the frequency dependence 
of the correlated noise. Since, the frequency variation should be rather weak, almost any 
smooth function with a only a few free parameters (polynomial, power law, gaussian, etc.) 
would probably work. 
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A. Leaky Modes 

The expansion of the wavefield into a discrete set of eigenmodes is appropriate only 
if the system lacks a continuous spectrum. Such spectra appear as branch cuts in the 
complex frequency plane. Here, I will consider two possible origins of such branch cuts. 
First, temporal damping of the form appearing in Equation (12. ip admits the possibility of 
overdamped solutions (OD). Such solutions manifest as a branch cut along the imaginary 
frequency axis. The existence of this branch cut is evident when one considers the frequency 
expression f2 = \fuP^^r^2A^. The square root function possesses branch points wherever 
its argument vanishes. Therefore, there are two branch points, one at w = and the 
other at a; = —2^7, with a corresponding branch cut connecting them along the imaginary 
axis. Frequencies along this branch cut correspond to overdamped waves whose temporal 
dependence is pure decay without oscillation. 

The second type of branch cut appears when the atmosphere possesses a global cut-off 
frequency, as the Sun does. If the product of the acoustic potential and the square of the 
sound speed remains bounded [c^V < a^) as the height becomes large {z — > oo), then waves 
of sufficiently large frequency {Q > a) are not trapped and leaky waves (LW) which radiate 
to infinity become possible. In the Sun, the condition for leaky waves is u)/27t > 5.2 mHz. 
There are two branch points associated with such waves, each located in the lower half of 
the complex frequency plane where fi^ = or w = ±-\/a^ — 7^ — i'y. I choose to connect 
these branch points seperately to infinity as shown in Figure HI and therefore the branch 
cut is bipartite. Frequencies along this branch cut correspond to propagating waves which 
radiate upwards through the photosphere, forming the high-frequency, pseudo-mode portion 
of the Sun's acoustic spectrum. 

When leaky waves are allowed by the boundary conditions, the system will not possess 
an infinite number of discrete eigenmodes. Instead, for a particular horizontal wavenumber 
k, the acoustic cavity will only have a finite number of discrete, trapped modes and a 
continuous spectrum of leaky modes. As before the trapped modes correspond to a zero in 
the Wronskian W, or equivalently to a simple pole in the Green's function. The leaky modes 
on the other hands correspond to the solutions along a branch cut. In an eigenfunction 
decomposition both types of waves are needed, as the discrete trapped waves on their own 
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do not form a complete set. For the problem at hand, these separate contributions can be 
calculated by Fourier transforms and contour integration in the complex frequency plane, 
which I perform in this appendix. 

I begin by taking the inverse Fourier transform of the Green's function as expressed in 
Equation I^M), 



G{z,z';k,t) 



^ I dtu G{z, z'; k, u) e"*^* , 



J_ f ^ Jl{r,r';k,u) 

2nJ'''' W{k,uj) 
c 



(Al) 



The integral is a contour integral through the complex frequency plane along the contour 
G which lies initially along the real frequency axis — see Figure Hl^a). Since all waves are 
damped, all of the poles and branch cuts appear in the lower half of the complex frequency 
plane. Therefore, for t < one may freely deform the contour upward to infinity and the 
integrand will vanish exponentially for large imaginary frequency, leading to the conclusion 
G = 0. However, for t > 0, in order for the integrand to vanish once again as the countour is 
deformed to infinity, the deformation must be downwards into the lower half plane. Since the 
lower half plane contains the poles and branch cuts, the contour must be deformed around 
them and the integral acquires contributions from the integration around each singularity. 
Using the residue theorem to evaluate the integration around each pole, one obtains the 
contribution from each trapped mode. 



G{z,z';Kt) = -^ "^"^^/(i^t J (1 - ^„o) ^(t) e^-^ 



n=-N 



Cbc 

where H{t) is the Heaviside step function and the prime on the Wronskian denotes differen- 
tiation with respect to the frequency argument, 

W'^"^. (A3) 
duj 



-17- 



Further, for future convenience, I have replaced the dummy variable of integration with u 
and the integration contour is disjoint and surrounds all of the branch cuts Cbc = C'od+C'lw- 
The summation over the index n is over a finite number (2A^) of poles, which includes two 
poles per mode, one with a positive real part of the frequency and the other with a negative 
real part, uj±n = ±Wn — ij- For clarity of notation, there is no pole associated with the 
index n = 0, and I have therefore included the factor with the Kronecker delta to exclude 
the n = term from the summation. 

Now that I have the Green's function in the form of an eigenfunction expansion, the 
result can be compared with Equation (13.21) by taking the temporal Fourier transform and 
including the effects of the source. 



G{z, z'\ k, u) 



ip{z, z'\ k, u) 



In the equations above I have used the following Fourier transform identity, 

oo 

[ dt H{t) e*('^-'^")* = — ^ . (A6) 

J OJ -OJn 

— oo 

One can easily verify that modulo the normalization of the eigenf unctions. Equation (]A5|) 
is equivalent to the modal expansion appearing in Section [HI Equation (13.50 . except for the 
appearance of the integral around the branch cuts (continuous spectra of leaky waves) and 
a summation over only a finite number of trapped modes instead of a countable infinity. 

Finally, note that due to the symmetries of the Wronskian W and the function \E', the 
integral around the branch cut associated with overdamped waves can be shown by direct 
evaluation to vanish. This is true for any problem in which the overdamped waves are 
confined to a cavity of finite extent. In truly infinite domains the precise cancellation around 
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the contour does not occur and the overdamped solutions may play a role. Therefore, the 
contour integral need only include the spectrum of leaky waves, Cbc = C*lw- 



A.l. Leaky Square- Well Potential 



To illustrate the mathematics developed previously in this appendix, I extend the rudi- 
mentary guitar-string model presented in the main body of this paper. Instead of requiring 
that all waves reflect from the cavity boundaries, I will examine a square-well potential that 
permits high-frequency waves to escape out the top of the cavity. Consider a one-dimensional 
acoustic cavity of length L with constant sound speed and a global cut-off frequency. The 
acoustic potential V{z) is identically zero within the cavity itself, and takes on a constant, 
nonzero value above the top of the cavity, z > L, 



= i 2, 2 -t T 

1^ a^/c^ if z > L . 

At the lower boundary of the cavity, z = 0, the solution must vanish, and at the top of the 
cavity the solution must match smoothly onto either an evanescent solution (for Q < a) or 
an upward travelling wave (for Q > a). If I assume that the driving layer lies within the 
cavity itself {z' < L), the homogeneous solutions ip^ and ip^ which obey the proper boundary 
conditions are as follows: 



^-{z;u) = sm{Kz) (A8) 

M^^ = ( f* - 'III '! ' * ^ (A9) 

[ exp[-K{z - L)] if z > L , 

where the constants K and k are, respectively, the wavenumber within the cavity and the 
evanescence rate for the region above the cavity. 



K{u,) = ^M, (AlO) 



Note, if ^2 > the evanescence rate k, becomes purely imaginary. For frequencies lying 
on the top of the branch cut, I choose the sign of the square root in Equation flAlip such 
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that the imaginary part of k is positive and the solution in the regime above the cavity, 
z > L, corresponds to an outward propagating wave and a radiation boundary condition 
is satisfied. Remember, the Fourier convention. Equation (I2.6p . is consistent with solutions 
with the time dependence ~ e"*"^*. 

The Wronskian of these two homogeneous solutions is of course sinusoidal. 



W{cj) = -[K cos(KL) + K sin(KL)] . (A12) 

If I assume that the observations are made above the cavity, z > L, as occurs for helioseismic 
observations, the Green's function has the solution 

KcosiKL) + K{KL) ■ ^^^^^ 

Figure shows the power spectrum and its components that result for this model when 
one assumes that the source is white (i.e., Sq is constant), 

P(: A _ I c |2 sm'(i^^0 exp [-2k{z - L)] 

The effect of the height of observation, exp [—2n{z — L)] (green curve), is a monotonically 
increasing function of frequency. The contribution from the Wronskian, |W|~^ (the red 
curve), behaves differently above and below the global cut-off frequency. Below the cut- 
off, a Lorentzian spike occurs at each mode frequency; above the cut-off, the Wronskian 
term oscillates and decays. The blue curve shows the driving efficiency, sin^(i^'2;'). It too 
is oscillatory, but with more rapid frequency variation than the Wronskian because z' < L. 
The zeros in the driving efficiency lead to the asymmetry in the total power profile (black 
curve) . 

In the alternate representation of a modal expansion, the Green's function and the power 
spectrum take on the form. 
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where the eigenfrequencies Un (or equivalently wavenumbers Kn and evanescence rates 
Kn are defined by the roots of the Wronskian, which provides a transcendental equation, 



W(w„) = - [Kn cos(K„L) + /€„ sin(K„L)] = 
fin 

A' = — 

»- r), — 1 



(A17) 
(A18) 

(A19) 



The wavenumber and evanescent length with the tilde, K and ^, are those corresponding 
to the integration frequency a), i.e., K = K{u) and k = k^u). Finally, the derivative of the 
Wronskian, evaluated at the mode frequencies, was used in these expressions. 
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Figure [5]^b) shows the resultant power spectrum (thick black curve) as well as the power 
associated just with the modes Pmodes (violet curve) and with the continuous spectrum Plw 
(aqua curve). 
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Similar to Figure [3l the thick orange curve shows the incoherent sum of the power in the 
modes. 
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Fig. 1. — Power spectrum as a function of frequency for a harmonic degree of I = 190 and an 
azimuthal order of m = 0. The spectrum was calculated from Dopplergram data obtained by the 
GONG network. Ten years of data were separated into overlapping three month subintervals and 
the individual spectra for all subintervals were averaged together. The three panels show different 
frequency ranges around the p4 mode, (a) The global context showing all of the ridges that exist 
below the acoustic cut-off frequency (fac ~ 5.2 mHz). (b) A narrower frequency range showing just 
the Pi ridge and its nearest neighbors, (c) A narrow frequency band that demonstrates that the p4 
ridge is a blending of the I = 190, n = 4 mode and all of the spatial leaks from modes of nearby 
harmonic degree. Differences in degree up to ±2 are clearly visible. The frequency profile of the 
blended line profiles are obviously skewed or asymmetric. 
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Fig. 2. — Power spectrum as a function of frequency for waves with / = 20 and m = 0. The same 
data as in Fig. [1] were used. For low values of harmonic degree (as illustrated here), each mode 
peak is well-separated from the spatial leaks of nearby modes. The leaks from modes degrees up 
to ±2 are shown here. Each of the frequency profiles (for the target mode as well as the leaks) is 
skewed and asymmetric. 
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Fig. 3. — The power spectrum for our simple "guitar-string" model. The source height is 
z' = 0.95 L and the observation height is at z = 0.99 L. The damping time is small compared 
to the sound-crossing time, 7-L/c = 0.1. The thick black curve shows the power spectrum of the 
solution as a function of frequency, and can be represented in either of two equivalent ways, equa- 
tion (j2.17p or (j3.12p . The thin colored curves show the contributions from three separate factors in 
equation (12.170 : the effect of the height of observation (green curve), the reciprocal of the square of 
the Wronskian (red curve), and the driving efficiency (blue curve). The total power is the product 
of these three factors. The Wronskian embodies the modal structure, and everywhere there is a 
resonant mode frequency the Wronskian vanishes, providing a peak with a Lorentzian profile. The 
driving efficiency possesses zeros at frequencies lieing between the modes. Asymmetry arises be- 
cause the zeros of the driving efficiency do not fall symmetrically between the mode frequencies. An 
alternative interpretation is that the modes represented in equation (j3.12p interfere with each other 
to produce the asymmetry. The thick orange curve demonstrates that the line profiles becomes 
symmetric if one neglects mode interference and assumes that the total power is the incoherent 
sum of the power in all modes. 
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Fig. 4. — Complex frequency space through which the inverse Fourier transform integral is per- 
formed, (a) The integration contour (red line labeled C) origininally lies along the real frequency 
axis, (b) For t > 0, the contour is deformed downwards and picks up contributions from each 
of the singularities. The diamonds mark the location of the poles of the Green's function, which 
correspond to the acoustic cavity's complex mode frequencies. The green crosses denote the branch 
points that frame the branch cut (green dashed line) associated with overdamped waves. The in- 
tegration contour around this branch cut is shown in red and labeled Cqd- The branch points and 
cut associated with leaky waves are indicated in blue, and the contour around this branch cut is 
labeled Clw- 
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Fig. 5. — The power spectrum for an acoustic cavity corresponding to a leaky square- well potential 
with a cavity depth of L. The source height is z' = 0.95 L, the global cut-off frequency is chosen 
to be Q = 18.0 c/L, and the damping parameter is the same as in Figure [3l 7 = 0.1 c/L. The 
observation height is at z = 1.1 L, just above the upper boundary of the cavity, (a) The thick black 
curve shows the power spectrum as a function of frequency obtained through equation ()A14p . The 
colored curves have the same meaning as in Figure [3l The vertical dashed line indicates the global 
cut-off frequency a. Asymmetry arises because the zeros of the driving efficiency (blue curve) do 
not fall symmetrically between the mode frequencies, (b) The power spectrum (thick black curve) 
obtained through an eigenfunction expansion, equation ()A16p . The violet curve shows the power 
associated with just the modes, -Pmodes) and the aqua curve shows the power due to the continuous 
spectrum of leaky waves, Plw- The thick orange curve presents the incoherent sum of the power in 
the modes, -Pine- The asymmetry apparent in the power is the result of interference between modes 
and between modes and the leaky waves. 



